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ABSTRACT 



Of 

o 



It is known that scattering of radiation by circumstellar dust can strongly change the line profiles in stellar spectra. 
This hampers the analysis of spectral lines originating in the emitting regions of heavily obscured young stars. To 
calculate the line profile of the scattered radiation, we suggest to use the approximation of remote scattering particles. 
This approximation assumes that the scattering dust grains are at a distance from the star that is much larger than the 
characteristic size of the emitting region. Using this method, we calculated the line profiles of several simple models. 
They show the Ha line profiles of Herbig AeBe stars in the presence and absence of motionless or moving dust. 
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1. Introduction 

Some amount of radiation scattered by circumstellar (CS) 
dust is present in the radiation of many young stars and 
is the main source of their linear polarization (Bastien & 
Landtsreet, 1979). In most cases, the contribution of this 
radiation is very small and does not exceed a few percent. 
Exceptions are three groups of young stellar objects whose 
scattered radiation can substantially exceed the direct stel- 
lar radiation. First, these are young stellar objects (YSO) 
embedded in opaque gas and dust envelopes. According to 
the current classification, these are objects of Class I (White 
et al. 2007). The extinction toward some of them exceeds 
20^ (see, e.g., Beck et al. 1991). Young stars surrounded 
by circumstellar disks seen edge-on belong to the second 
group. In these stars, the direct radiation is strongly at- 
tenuated due to the extinction in the disk, and one can 
observe them only via scattered radiation (Padgett et al. 
1999). The T Tauri star HH 30 is a prototype of these 
objects (Burrows et al. 1996). The third group consists of 
the UX Ori type stars. The photometric activity of these 
young stars is caused by the nonhomogeneous structure of 
their circumstellar disks and their small inclination to the 
line of sight (Grinin et al. 1991; Natta et al. 1999). Opaque 
fragments of the CS disk intersect the line of sight from 
time to time. During these events, the direct radiation of 
the star decreases while the contribution of the scattered 
radiation increases, which is testified by the growth of the 
linear polarization of the star (Grinin et al. 1991). 

In all these cases, an observer obtains the information 
on the stellar spectrum by observing the radiation scat- 



Send offprint requests to: V. P. Grinin 



tered by the circumstellar dust. It is well known that in 
the coordinate system of the dust particle, the photon fre- 
quency before and after scattering is the same. However, 
the resulting scattered-light spectrum is not the same as 
the spectrum in the absence of scattering. 

As noted by Appenzeller, Bertout & Stahl (2005), scat- 
tered light is not a precise source of information about 
obscured objects. These authors observed several T Tauri 
stars with edge-on disks. The analysis of their spectra 
showed that the disks are completely opaque at visible 
wavelengths and that light from the central objects reaches 
an observer only via scattering layers above and below the 
disk planes. According to these authors, light could be scat- 
tered either by matter located along the rotation axis or in 
the disk atmosphere. As a result, the spectrum of scattered 
radiation can differ from the spectrum of the CS emitting 
region. Grinin, Mitskevich & Tambovtseva (2006) consid- 
ered the broadening of photospheric lines due to the scat- 
tering by the CS dust at the puffed- up inner rim of the 
accretion disks. This effect can be observed in UX Ori type 
stars during eclipses by opaque fragments of CS disks. The 
transformation of the spectral line profiles due to dust scat- 
tering in expanding dusty envelopes can also be observed 
in the spectra of evolved stars (see Lefevre 1992, and refer- 
ences therein). 

Below we consider several model situations typical for 
embedded YSOs (Fig.[T]) and show how the profiles of spec- 
tral lines broadened by motion in the emitting gas (e.g., disk 
wind) can change after scattering by circumstellar dust par- 
ticles. 
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Fig. 1. Sketch of a deeply embedded young star with a ro- 
tating accreting disk, disk wind, and an envelope (cocoon) 
with low-density polar cavities. Dust particles are assumed 
to be near the polar axis and walls of the cavity. Possible 
profiles of emission lines are shown for several viewing con- 
ditions: direct (A) and scattered (B and C) radiation from 
the central object, which consists of the star, and its sur- 
rounding emitting region (e.g., disk and disk wind region). 
Here /g is the emergent radiation from this central object, 
and is the radiation scattered by the dust particles 
in cases B and C respectively, r^, tb and rc are the optical 
thickness of the circumstellar envelope at different heights 
above the disk {tb <C ta and rc ^ t^), ^ is the viewing 
angle, i' is the angle between the polar axis and direction 
to the scattering particle, and 6> is a geometrical parameter. 



2. Scattering by the motionless dust 

The effect of scattering by dust particles depends on the ge- 
ometry and kinematics of the line-formation region and the 
scattering region (see Fig. [T]). If the spectral line arises, for 
example, in a flattened rotating envelope or accretion disk, 
the line profiles of the direct and scattered radiation can 
be very different: observer A, for example, would observe 
a wide double-peaked or P Cygni line profile from the ac- 
cretion disk and/or disk wind. However, if the direct view 
to the emitting region is blocked, then observer B would 
observe a narrow single or double-peaked line. If there 
scattering dust near the polar axis of the disk far aw 
from the disk, then observer C would observe a very m 
row line, since the radial motion seen at the position of t 
scattering dust is almost zero. These three simple examp] 
show that the observed spectrum depends on both the ^ 
ometry and kinematics of the emitting region itself and t 
distribution of the scattering particles in the neighborho* 
of the star. 

In the general case when the spectral line originates ir 
moving gas with an arbitrary velocity field, the expressL 
for the intensity /sc(^) at the frequency u of the radiati< 
scattered at point A (Fig. [2]) by a single motionless (in t 
coordinate system of the star) particle is 

Isc{i^)=7ra^QscJ 5(r)(/)(r,^-^o^) X ( 
V 



Here V is the volume occupied by the emitting gas, a is 
the radius of the particle, Qsc, and 0' are the efficiency 
factor for scattering, phase function, and the scattering an- 
gle, respectively, S{r) is the source function for the transi- 
tion considered at the point r of the emitting region, (j) is 
the profile of the absorption coefficient normalized to the 
unity, Vi^(r) is the projection of the velocity of the emitting 
gas in the point r onto the vector R, which connects this 
point with the scattering particle at point A, 

rd = r^ + T^', (2) 

where is the optical thickness between point r and point 
A caused by the CS extinction, r'J is the optical thickness 
between point A and the observer; hz is the integrated line 
opacity in the considered spectral line, and r{u^ r) is the 
line optical depth at point r and for the frequency u in the 
direction of point A 

r(^,r)= ,^{r')^{jy - jy^^^^i^)dR' . (3) 
Jo c 

Here the integration is performed along the vector R linking 
the scattering particles with the point r; r' and are the 
integration variables. We assume a complete redistribution 
in the line frequencies in the reference frame of the atom 
and use the Doppler profile (j) in the calculations of r and 
the line intensity. 

In order to take into account scattering by all parti- 
cles, we have to integrate the expression (1) over the to- 
tal volume occupied by the dust, as well as to consider all 
sizes and types of dust particles participating in scattering. 
The solution of this problem faces substantial difficulties. 
Therefore, below we consider the simplified model when 
scattering grains are at some fixed distance D = |D| from 
the star, and D is much larger than the characteristic size of 
the emitting region. This approximation is quite realistic in 
many cases since the region of the permitted line formation 
is usually much more compact than the region occupied by 
the CS dust. 

2.1. Approximation of remote scattering particles 

In this case, we can use the following simplifications: i) in 
integral (1), it is possible to replace hy and thus 




Fig. 2. Geometric parameters of the considered model (a 
central cross-section of the CS region). 
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to extract from the integral; ii) in the phase function 
/(6>'), it is possible to replace 6' by V^, where i/j is the angle 
between vector D and direction from the point A to the 
observer, and also extract it from the integral. We can do 
the same with the optical thickness Td if we use the optical 
thickness between the points O (location of the star) and 
A instead of r^. After these simplifications we obtain the 
asymptotic equation for single scattering by the elementary 
volume at point A, which can be written in the following 
form: 



(4) 



where N is the number density of particles, /oo(^) is the over 
the whole emitting region integrated intensity at frequency 
u in the direction to the remote observer whose line of sight 
coincides with the vector D: 



(5) 



Since the scattering particle is assumed to be motionless in 
the star coordinate system, the frequency of the scattered 
radiation in any direction will coincide with the frequency 
of the incident radiation, and since the values of /{ip) and 
Qsc before the integral are practically constant within the 
profile of the spectral line, the value /oo(^) will determine 
the line profile in the spectrum of the scattered radiation. 
It is important that this profile is the same for all scattering 
angles and depends only on the angle . 

It should be stressed that the intensity of the radiation 
in the expression given above corresponds to the case of sin- 
gle scattering. Since we assume here that the dust particles 
are motionless in the star coordinate system, the photon 
frequency will be conserved in the next scattering. Thus, 
the first scattering plays a key role in the transformation of 
the line profile in the considered case. It should be noted, 
that in the optical region of the spectrum, albedo of scat- 
tering by CS dust is about 0.5 (see, e.g., Natta & Whitney 
2000). Therefore, the first scattering gives the main contri- 
bution to the intensity of the scattered radiation. 
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Fig. 3. The profile of the Ha line in the spectrum of the 
radiation scattered at point A (a,b,c,d) and the spectrum 
of direct radiation (e). The parameters of the wind model 
are described in the text and the Appendix. The angle 
between the symmetry axis of the disk and the vector D are 
5° (a), 20° (b) ,40° (c), and 60° (d). The angle between the 
line of sight and the symmetry axis of the disk is i = 80° 
(e). 

shows the Ha line profile which the observer could see in 
the absence of scattering. We see that this profile (broad- 
ened by the motion of the emitting gas) has a P Cygni 
shape (instead of pure emission) and it is wider than the 
profiles observed after scattering at the point A. This ex- 
ample shows that even in the case of the motionless dust, 
the profile of the spectral line of the direct radiation and 
the radiation after scattering can strongly differ. 



2.2. The Ha line in the spectrum of the scattered radiation 

As an example, let us consider the following model: the 
emission spectrum of the young star is formed in the axi- 
ally symmetric disk wind, which starts from the surface of 
the accretion disk. The star and the emission region are sur- 
rounded by an optically thick axially symmetric envelope 
with cavity (see Fig. 1), whose polar region is transparent 
for the radiation. The spectrum of the radiation scattered 
by the dust on the surface of the cavity walls and by dust 
in the polar region is observed as the spectrum of the star. 

For the calculations of ifci^) of the Ha line, we used 
the model of the magnetocentrifugal disk wind of Ae/Be 
Herbig stars reported by Grinin & Tambovtseva (2011). A 
description of this model is given in the cited paper and the 
Appendix. 

The profiles of the Ha line in the spectrum of the scat- 
tered radiation are calculated for several cases where the 
angle between the symmetry axis of the wind and the 
direction to the scattering particle ranged from 5 to 60 de- 
grees (Fig. [3t^-d), and the angle between the line of sight 
and the symmetry axis of the disk was i = 80°. Figure [3^ 



3. Scattering by the moving dust 

In the reality, the dust particles are moving. The motion 
may be purely Keplerian rotation of the CS disk or purely 
radial (bi-conical) outfiow, or their combination, like in the 
dusty disk wind (Safier 1993; Konigl & Salmeron, 2011). 
Scattering of photons by moving dust provides an addi- 
tional Doppler shift of the photon frequency, which leads 
to an additional deformation of the line profiles. In this 
case, in the equation dU we have to replace u hy ^ where 

= v -\- [vs -\- Us — V COS {6 — i')]/c (6) 

Here Vg and Us are the projections of the radial {v) and 
rotation (u) velocity components at point A on the line of 
sight s: 

Vs = v{cosO cosi -\- sinO sini coscp) ^ (7) 

Us = —usinismcp^ (8) 

where cp is the azimuthal coordinate angle of point A (see 
Fig.Sl). 
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The term —vcos {0 — i')/c in Eq. (j6j) takes into account 
the Doppler shift of the frequency v due to the motion of 
the scattering particles relatively to the emitting region. 



Obs. 




Fig. 4. The sketch of a dust particle at point A moving 
with the velocity components v and u and scattering the 
photons in the direction to an observer. 



The intensity of the scattered radiation at the frequency 
V in the direction to an observer is the integral over the 
whole volume V: 

hM= I lt{v')dY (9) 
Jv 

Below we consider two limiting cases: a) radial motion of 
the dust particles with constant velocity v and tangential 
velocity = 0, and b) tangential motion of the dust with 
the velocity u and radial velocity v = 0. The intensity of 
the scattered radiation from all elementary segments of the 
scattering volume with fixed coordinates 6>, D is: 

/,,(^)(x r ii{v')d^, (10) 

where /^(i^) is determined from equation (jH but with 
changes described above (Eqs.® - (|8|)). 

When calculating the line profile of the scattered radia- 
tion on the base of this formula, we assumed for simplicity 
that = i\ the optical thickness between the scattering 
particles and the observer r" <C 1, and the scattering is 
isotropic. 

The results of the calculations for the limiting cases de- 
scribed above are presented in Figs. [5] and [6l According to 
Fig. O a radial motion of the dust from the star provides 
a red shift of the line profile. This shift will be larger if 
the main contribution to the scattered light is provided by 
the back wall of the conical cavity (Fig. [6]), which is ob- 
served in some heavily obscured stars (see, e.g., Preibisch 
et al. 2003). The rotation does not shift the line profile, but 
smooths it and makes it broader (FiglSb). 

These examples show that the effect of motion of the CS 
dust leads to more divers transformations of the line profiles 
compared to the case of motionless dust. In the both cases 
considered above, the line profile of the scattered radiation 
strongly deviates from the profile of the direct radiation in 
the absence of scattering presented in FiglSt. 
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Fig. 6. The profile of the Ha line in the spectrum of the 
scattered radiation with (solid) and without (dots) radial 
outfiow of the dust with v = 100 kms~^. The contribu- 
tion of the back wall is only considered here. Therefore, the 
angle cp in Eq. ([8]) is varied from 7r/2 to 37r/2. The other 
parameters are the same as in Figure [5l See text for more 
details. 



4. Summary and concluding remarks 

The presented calculations show that the profiles of spectral 
lines arising in moving CS gas can change when scattering 
by CS dust takes place. The distortion can appear even if 
the CS dust is not moving in the coordinate system of the 
star. The lines have more diverse forms in the case of the 
moving dust. This is important to keep in mind when an- 
alyzing emission spectra of heavily obscured young stars, 
and it is especially important for the comparison of spec- 
tral lines at different wavelength (e.g., Ra and Br7), since 
the dust extinction depends on the wavelength. In heav- 
ily obscured young objects, the Ha line is observed via the 
scattering radiation in practically all cases. The contribu- 
tion of the scattered light to the Br7 line is less important 
if Ay < 7 — 8^. However, in young objects with Ay ^ 8^, 
even the Br7 is predominatingly observed in the scattered 
light and its profile can also be distorted by scattering. In 
these cases, only observations of mid- and far-infrared spec- 
tral lines may avoid these scattering effects. 

The changes of the spectral fine profiles due to the scat- 
tering by the CS dust is also valid for the photospheric 
spectrum of the star itself since its radiation propagates in 
a similar same way to an observer as the radiation of the 
CS emitting region. As a result, the measured radial and 
rotational velocities of the heavily obscured star can differ 
from their real values. 

It should be noted, that in the case when the spectral 
line is formed in a spherically symmetric stellar wind, the 
spectrum of radiation scattered by a single particle at an 
arbitrary remote point will be the same for any scattering 
angle and will coincide with the spectrum of the direct ra- 
diation of the star. This is the only kind of model where 
the profiles of the spectral lines do not change in the case 
of dust scattering. 

Finally, we emphasize that highly obscured YSOs are, 
of course, not the only objects in which the scattered ra- 
diation can be an important component of the observed 
radiation. Similar conditions exist in active galaxy nuclei 
observed under certain inclinations (see Antonucci & Miller 
1985, Antonucci 1993; Cohen et al. 1999; Smith et al. 2005, 
and references there). In these objects, the scattered radia- 
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Fig. 5. The profile of the Ha hne in the spectrum of the scattered radiation with (solid) and without (dots) assuming dust 
motion: a) line profile transformation due to a radial outflow with velocity v = 100 kms~^; b) line profile transformation 
due to rotation {u = 100 kms~^); c)line profile of direct radiation in the absence of dust scattering. The angle between 
the symmetry axis and the direction to the scattering particle is = 40°, the angle between the symmetry axis of the 
disk and the line of sight is i = 80°. 



tion arises due to the Thomson scattering on free electrons 
and also due to scattering by dust particles in the dusty 
torus and the Narrow Line region. Spectropolarimetric ob- 
servations have played an important role in studies of the 
physical properties of these objects. Similar spectropolari- 
metric observations would also be very useful in studies of 
heavily obscured YSOs. 

Acknowledgements. This research has been supported in part by the 
program of the Presidium of RAS N 21 and grant N. Sh.- 1625.2012.2. 



References 

Antonucci, R.R.J. 1993, ARA&A 31, 473 
Antonucci, R.R.J, & Miller, J.S. 1985, ApJ 297, 621 
Appenzeller, I., Bertout, C, & Stahl, O. 2005, A&A 434, 1005 
Blandford, R.D., &; Payne, D.G. 1982, MNRAS 199, 883 
Bastien, P. &; Landstreet, J.D. 1979, ApJ 229, L137 
Beck,S.C., Fischer, J., Smith, H.A. 1991, ApJ 383, 336 
Burrows, C.J., Stapelfeldt, K.R., Watson, A.M., et al. 1996, ApJ 
473, 437 

Cohen, M., Ogle, P.M., Tran, H.D., et al. 1999, AJ 118, 1963 
Grinin, V.P., Kiselev, N. N., Chernova, G. P. et al. 1991, Ap. Sp. 
Sci. 186, 283 

Grinin, V.P., Mitskevich, A.S., Tambovtseva, L.V., 2006, Astronomy 
Letters 32, 110 

Grinin, V.P., & Tambovtseva, L.V. 2011, Astronomy Reports 55, 
704 

Konigl, A. &; Salmeron, R. in Physical Processes in Circumstellar 
Disks around Young Stars, Ed. Paulo J.V. Garcia. Chicago, IL: 
University of Chicago Press 2011, p. 283 

Kurosawa, R. Harries, T. J., & Symington, N. H. 2006, MNRAS, 
370, 580 

Lefevre, J. 1992, A&A 254, 274 

Malbet, F., Benisty, M., de Wit, W.-J. et al. 2007, A&A 464, 43 

Natta, A., Whitney, B., 2000, A&A, 364, 633 

Natta, A., Prusti, T., Neri, R. et al. 1999, A&A 350, 541 

Padgett, D.L., Brandner, W., Stapelfeldt, K. R. et al. 1999, AJ 117, 

1490 

Preibisch, T., Balega, Y.Y., Schertl D., & Weigelt, G. 2003, A&A 
412, 735 

Safier, P. 1993, ApJ 408, 115 

Smith, J.E., Robinson, A., Young, S., et al. 2005, MNRAS 359, 846 
White, R.J., Greene, T.P., Doppmann, G.W. et al. 2007, in 
Protostars & Planets V., B. Reipurth, D. Jewitt, and K.Keil (eds), 
University of Arizona Press, Tucson, 2007, p. 951 
Weigelt, G.,Grinin, V. P.; Groh, J. H.; et al. 2011, A&A 527, 103 



5 



Grinin et al. 



Appendix: the disk wind model 

For the description of the disk wind, we use a coordi- 
nate system 6>, uj) centered at point S (see the sketch in 
Fig. [7j). In MHD models of the disk wind, the inchnation 
angle of the first streamline with respect to the disk plane 
is typically assumed to be around 60° (Blandford & Payne 
1982). However, in the case of the hot Herbig Be stars, the 
strong stellar radiation pressure can bend down the stream- 
lines and make the disk wind flatter. The kinematical model 
described below is typical for disk- wind models adopted for 
different astrophysical objects with accretion disks includ- 
ing young stars stars (see e.g. Kurosava et al. 2006), and 
references there). It is assumed that the tangential velocity 
of the wind changes along streamlines as: 

u{uj) = uxiuJi) {uj/oJi)~^ (11) 

where cu = I sinO is the distance of the point (/, 0) from the 
rotation axis and uji = U sin6>; uk {(^i) = {G M^/uji)^^'^ at 
the point {uJi) at the base of streamline i . 

The radial velocity v increases along the streamlines as 

v{l) = ^0 + (^oo - ^o) (1 - h/lf (12) 

where '^o and ^oo are the initial and terminal values of the 
radial velocity and /3 is a parameter. We adopt for simplicity 
'^oo = f uxioji)^ where uxioJi) is the Keplerian velocity at 
distance cui from the disk axis, and / is a parameter. As in 
the Kurosawa et al. (2006) paper we used / = 2. 

The local mass- loss rate per unit area of the disk, m^^, 
is a function of the cylindric radius uj. For the description 
of m^y, we use a simple power law 

m^{uj) ^ uj~^ , (13) 

where 7 is the input parameter. The other input parameter 
is the total mass loss rate integrated over all the disk wind 
ejection region (from both sides of the disk) M^. 

Using the continuity equation, we calculated the dis- 
tribution of the number density in each point of the disk 
wind. The calculations of the ionization state and the num- 
ber densities of the atomic levels were performed for the 
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Fig. 7. Sketch of the geometry of the disk-wind model 
adopted. 



isothermal disk wind model with the account of both col- 
lision and radiative processes on the base of the Sobolev 
approximation (see Grinin & Tambovtseva (2011) for more 
details). 

The radiation of the disk wind was calculated for the 
following model parameters: the mass loss rate = 
10-^ Mq yr-i; the electron temperature Te = 8000 K, 

= 2.5M0, = 2.4i?0, 7 = 1, the disk wind launching 
region is located between 0.1 and 1 AU. The half opening 
angle of the disk wind 61 = 60° (see Fig. [7j). The radia- 
tion of the star is described by the Kurucz model with the 
effective temperature Tg/ = 10^ K and log^ = 4. 
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